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Abstract 

Using the constrained curve fitting method and overlap fermions with the 
lowest pion mass at 180 MeV, we observe that the masses of the first posi- 
tive and negative parity excited states of the nucleon tend to cross over as the 
quark masses are taken to the chiral limit. Both results at the physical pion 
mass agree with the experimental values of the Roper resonance (A^ 1 / 2+ (1440)) 
and S*n (A fl / 2_ (1535)). This is seen for the first time in a lattice QCD calcu- 
lation. These results are obtained on a quenched Iwasaki 16 3 x 28 lattice with 
a = 0.2 fm. We also extract the ghost n'N states (a quenched artifact) which are 
shown to decouple from the nucleon interpolation field above m n ~ 300 MeV. 
From the quark mass dependence of these states in the chiral region, we con- 
clude that spontaneously broken chiral symmetry dictates the dynamics of light 
quarks in the nucleon. 



The Roper resonance has been studied extensively, but its status as the first 
excited state of the nucleon with the same quantum numbers is intriguing. Although 
awarded four stars in the Particle Data Table, it is not as easily identifiable as other 
resonances. It does not show up as a peak in the cross sections in various experiments; 
rather its existence is revealed through phase shift analysis 0121- Furthermore, the 
theoretical situation is in a quandary. First of all, it has been noted for a long time 
that it is rather unusual to have the first positive parity excited state lower than the 
negative parity excited state which is the A 1//2 ~(1535) in the Sn irN channel. This 
is contrary to the excitation pattern in the meson sectors with either light or heavy 
quarks. This parity reversal has caused problems for the otherwise successful quark 
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models based on £77(6) symmetry with color-spin interaction between the quarks |3] 
which cannot accommodate such a pattern. Realistic potential calculations with 
linear and Coulomb potentials |3] and the relativistic quark model jH] all predict 
the Roper to be ~ 100 - 200 MeV above the experimental value with the negative 
parity state lying lower. On the other hand, the pattern of parity reversal was readily 
obtained in the chiral soliton model like the Skyrme model via the small oscillation 
approximation to irN scattering [Hj. Although the first calculation [7j of the original 
skyrmion gives rise to a breathing mode which is ~ 200 MeV lower than the Roper 
resonance, it was shown later [Hj that the introduction of the sixth order term, which 
is the zero range approximation for the u meson coupling, changes the compression 
modulus, yielding better agreement with experiment for both the mass and width in 
ttN scattering. 

Since the quark potential model is based on the 577(6) symmetry with residual 
color-spin interaction between the quarks, whereas the chiral soliton model is based 
on spontaneous broken chiral symmetry, their distinctly different predictions on the 
ordering of the positive and negative parity excited states may be a reflection of dif- 
ferent dynamics as a direct consequence of the respective symmetry. This possibility 
has prompted the suggestion [H] that the parity reversal in the excited nucleon and 
A, in contrast to that in the excited A spectrum, is an indication that the inter-quark 
interaction of the light quarks is mainly of the flavor-spin nature rather than the 
color-spin nature (e.g. one-gluon exchange type.) This suggestion is supported in the 
lattice QCD study of "Valence QCD" [THj where one finds that the hyperfine split- 
ting between the nucleon and A is largely diminished when the Z-graphs in the quark 
propagator are eliminated. This is an indication that the color- magnetic interaction 
is not the primary source of the inter-quark spin-spin interaction for light quarks. 
(The color-magnetic part, being spatial in origin, is unaffected by the truncation of 
Z-graphs in Valence QCD, which only affects the time part.) Yet, it is consistent 
with the Goldstone-boson-exchange picture which requires Z-graphs and leads to a 
flavor-spin interaction. 

The failure of the S77(6) quark model to delineate the Roper and its photo- 
production has prompted the speculation that the Roper resonance may be a hy- 
brid state with excited glue [TT] or a qqqqq five quark state [12]. Thus, unraveling 
the nature of Roper resonance has direct bearing on our understanding of the quark 
structure and chiral dynamics of baryons, which is one of the primary missions at 
experimental facilities like Jefferson Lab. 

Lattice QCD is perhaps the most desirable tool to adjudicate the theoretical con- 
troversy surrounding the issue. In fact, there have been several calculations to study 
the positive-parity excitation of the nucleon [T3J El EES UH UZ1 UE] • However, they 
have not been able to probe the relevant low quark mass region while preserving chi- 
ral symmetry at finite lattice spacing (except Ref. ^3] which uses the Domain Wall 
fermion). We employ the overlap fermion JH] on a large lattice which admits cal- 



2 



CO 

.2 

o 

Pi 
.2 

o 

S-H 
S-H 

O 

O 



CO 



3e-06 
2e-06 
le-06 


-le-06 
3e-06 

2e-06 

le-06 



-le-06 



- 

: 1 


= 0.342(6) 


;1 


= 0.290(6) 




m, = 0.248(7) 


i 








i 


J 

- 


= 0.212(7) 

i 


; 1 


m» = 0.196(7) 


= 1 


rn v = 0.181(8) 





i i 


\^ 




i i 


a— 

-- 



2 4 6 8 10 12 14 2 4 6 8 10 12 14 2 4 6 



10 12 14 16 



Time-Slices (t) 



Figure 1: Sn correlators for six low quark masses. Negative dip of the correlators is 
an indication of the domination of the ghost rj'N state over the physical Sn. are 
in GeV. 



culations with realistically small quark masses [20] to study the chiral region. Since 
the controversy about the nature of Roper hinges on chiral symmetry, it is essential 
to have a fermion action which explicitly exhibits the correct spontaneously broken 
chiral symmetry. 

Another difficulty of the calculation of the excited states in lattice QCD is that the 
conventional two-exponential fits are not reliable. Facing the uncertainty of the fitting 
procedure for the excited state, it has been suggested to use a non-standard nucleon 
interpolating field which vanishes in the non-relativistic limit El EE EH CBj , in 
the hope that it may have negligible overlap with the nucleon so that the Roper state 
can be seen more readily. However, the lowest state calculated with this interpolation 
field (2.2 GeV) [TIIlEllUnillEllIE] is much higher than the Roper state. Employing the 
maximum entropy method allows one to study the nucleon and its radial excitation 
with the standard nucleon interpolation field |17j . However, with the pion mass at 
~ 600 M eV, the nucleon radial excitation is still too high (~ 2 GeV) . In order to 
have a reliable prediction of the excited state in lattice QCD calculations by way of 
addressing the above mentioned difficulties, we implement constrained curve fitting 
based on an empirical Bayes method. It has been advocated [2H 122] recently as a 
powerful tool which utilizes more data points while better controlling the systematics. 
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Since we adopt the quenched approximation, we should mention that there are 
quenched artifacts. Due to the absence of quark loops in quenched QCD, the would-be 
rj' propagator involves only double rj poles in hairpin diagrams. This leads to quenched 
chiral logs in hadron masses which are clearly observed and extracted from our recent 
lattice calculation of pion and nucleon masses j2Hj- Another quenched artifact is the 
contribution of ghost states in the hadron correlators. As first observed in the a® 
meson channel [%4*j . the ghost S- wave rj'n state lies lower than ao for small quark 
mass. There is a similar situation with the excited state of the nucleon where the 
P-wave rj'N appears in the vicinity of the Roper. Since this is not clearly exhibited 
in the nucleon correlator where the nucleon is the lowest state in the channel and 
dominates the long-time behavior of the correlator, we can look at the parity partner 
of the nucleon (Nz~~ or Sn) with / = 1/2. There, the lowest S-wave rj'N state 
with a mass close to the sum of the pion and nucleon masses can be lower than Sn 
for sufficiently low quark mass. Since the 77 - rj coupling in the hairpin diagram is 
negative j21], one expects that the Sn correlator will turn negative at larger time 
separations as in the case of the a [211 • I n Fig. E we show the Sn correlators for 
6 low quark cases with pion mass from m n = 181(8) MeV to m n = 342(6) MeV. We 
see that for pion mass lower than 248(7) MeV, the Sn correlator starts to develop a 
negative dip at time slices beyond 4, and it is progressively more negative for smaller 
quark masses. This is a clear indication that the ghost rj'N state in the S-wave is 
dominating the correlator over the physical Sn which lies higher in mass. 

One can calculate the one-loop hairpin diagram on a L 3 x T lattice which gives 
the following contribution to the zero-momentum nucleon correlator for the P-wave 
would-be rj' loop 

.„ /jA 2 \Z\ 2 T ^ E N + E w + m N 

AG NN (t) = -m 0W E 

p 2 



m 



[1 + E n t)e- {EN+E ^. 



N 



Here — is the coupling between the would-be Goldstone bosons. However, this is 
not quite adequate to be used in the fitting because there is an interaction between 
the would-be rj' and the nucleon which can lead to an energy which is different from 
En + E n . Instead of modeling it with the re-summation of the hairpin diagram, as 
is done for the ao meson (21], we shall introduce these ghost states which are in the 
vicinity of the Roper and Sn with the form W(l + E v t)e~ E ^' Nt . This has the same 
functional form as in Eq. The (1 + E n t) pre-factor is preserved to reflect the 

double-pole nature of the hairpin diagram and the weight W is constrained to be 
negative to preserve the ghost nature of the state. The only exception is E v >n which 
is allowed to be different from but close to the non- interacting energy E^ + E^. Being 
the energy of the interacting would-be rj' and N state, it will be fitted to the data 
together with the negative weight W. Since we work in a finite box, the rj'N states are 
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Figure 2: Roper (left) and Su (right) masses with their associated rj'N ghost states. 



discrete and they will be constrained to be near the energy of the two non-interacting 
particles, with discrete lattice momenta p n = =7p, n — 0, ±1, ■ ■ -. 

With a constrained-curve fitting method j22| . we develop an adaptive sequential 
way of extracting the priors from data at a smaller volume, i.e. 12 3 x 28. We fit the 
ground state first over a few large-time slices, and then use the fitted values as priors 
to constrain the ground state in a subsequent fit (which includes earlier time slices) 
of the excited state. We include the next excited states, one by one, constraining all 
previously fitted states to their most recently fitted values. For each new fit, we add 
one more earlier time slice and check if it is needed by comparing its contribution to 
the statistical error of the added time slice. When the time slices are exhausted, we 
use the resultant fit parameters as priors for a standard constrained-curve fit of the 
correlator on the larger 16 3 x 28 lattice. More details can be found in Ref. [22] • For 
the rj'N state, we make an appropriate adjustment to priors for E^n and E n due to 
the fact that the discrete momentum is inversely proportional to the spatial lattice 
size L. We find that the constraint to make W negative is crucial to obtain a fit with 
acceptable x 2 ; especially for the case of the nucleon correlator where the Roper and 
the first rj'N state are nearly degenerate. 

We shall report our results on 16 3 x 28 and 12 3 x 28 lattices with quenched Iwasaki 
gauge action and overlap fermions with a lattice spacing a = 0.200(3) fm determined 
from f 7T (m 7T ). This makes our lattice sizes as large as 3.2 fm and 2.4 fm. At our lowest 
pion mass at 181(8) MeV, the finite volume error is estimated to be ~ 2.7% We 
first present the Roper and Su masses for the 16 3 x 28 lattice with their associated 
ghost rj'N states in the P-wave and S'-wave respectively They are plotted as a 
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Figure 3: Ratio of spectral weights from 12 3 to 16 3 lattices for Roper and P-wave 
rj'N (left) and Sn and S-wave rj'N (right). 



function of in Fig. |2l We see in Fig. |2] (left) that the rj'N state (p = 2n/L) and 
the Roper are nearly degenerate; they are resolved when the weight for the rj'N is 
constrained to be negative. Fig. El (right) shows that the p = S-w&ve rj'N lies lower 
than Sn, and below m n ~ 316 MeV the p = 2n/L rj'N state also becomes lower than 
Sn- The dashed line in Fig. |2] (left) is the mass for the non- interacting would-be rj' 
and N which lies below the fit of the rj'N state. This shows that there is a repulsive 
interaction of about 50 MeV between the would-be rj' and iV in the P-wave. The 
corresponding non-interacting rj'N states, represented by the solid line (p = 0) and 
the dashed line (p = 2n/L) in Fig. |21 (right), are closer to the interacting rj'N states 
with ~ 10 — 50 MeV repulsion in this case. 

To verify that the rj'N state can be separated from the physical Roper and Sn, we 
examine the ratio of the spectral weights for the Roper (respectively, the Sn) and the 
rj'N from the 12 3 x 28 lattice to that on the 16 3 x 28 lattice. (See Fig. El) The ratio 
for the Roper (Sn) is consistent with unity which is expected for a one-particle state; 
whereas, the corresponding ratios for the rj'N states are more consistent with the ratio 
of the three volumes Vs(16) /Vs(12) = 2.37 which reflects the volume dependence of 
the two-particle scattering state. 

In Fig. HI we display the weight of the rj'N state and those of the Roper and 
Sn. The weight for the Roper is essentially constant (Fig. HI (left)) and that of 
Sn decreases as the pion mass increases and reaches a constant for m n > 400 MeV 
(Fig. HI (right)). The negative weight of the rj'N is 3 orders of magnitude smaller than 
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Figure 4: Weights of Roper and rj'N (left), and S\i and rj'N (right). They are in 
arbitrary units. 



the Roper above m n = 266 MeV in the nucleon channel and = 320 MeV in the 
Su channel. This decoupling means that the two-particle intermediate state rj'N is 
not produced by the three-quark nucleon interpolating field xjy so that the matrix 
element (?/^V|xjv|0) 1S zero - This confirms our earlier study of Valence QCD JU] which 
suggests that, for higher quark mass, the SU(6) symmetry prevails which explains 
why the valence quark model with color-spin interaction works well. However, for 
light quarks, it is chiral symmetry which dictates the dynamics via the Z-graphs. 
This is to be manifested as the meson loops as well as meson exchanges between 
quarks in the chiral quark model of the baryon [23] . The present study suggests that 
this transition occurs in the region m n ~ 300 MeV. We should point out that the fact 
that the rj'N ghost states are separated out from the physical Sn(1535) in our fitting 
relies on the fact that the ghost state has a different functional form from that of the 
physical state and its weight is negative. When the dynamical fermion calculation is 
available with realistically light quarks, these rj'N ghost states will turn into physical 
rjN and rj'N states. Since Sn (1535) is close to the rjN threshold and both have the 
same functional form, i.e. W e~ mt , it would be difficult to separate them with a fitting 
procedure based on one interpolation field. 

Finally, we plot the nucleon, Roper and Sn in Fig. El as a function of m%. We 
see that for heavy quarks (m n > 800 MeV), the Roper, Sn, and nucleon splittings 
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Figure 5: Nucleon, Roper, and Su masses in GeV as a function of m 2 . The experi- 
mental values are indicated by the corresponding open symbols. 

are like those of the heavy quarkonium. When the quark mass becomes lighter, the 
Roper and Su have a tendency to coincide and possibly cross over. More statistics 
are needed to clarify this point. However, we have plotted as the insert in Fig. the 
ratio of Roper to nucleon as a function of m^(GeV 2 ) which has a smaller error (by 
~ 40%) than the Roper mass itself. It shows that the Roper is consistent with the 
experimental value near the physical pion mass. The chiral behavior of the nucleon 
mass is analyzed in more detail with the inclusion of the m 7r log(m 7r ) term in |23j . 
Here, we shall simply use the form C + Ci/2m n + Cim^ to extrapolate the Roper and 
Su to the physical pion mass. The resultant nucleon mass of 928(56) MeV, Roper at 
1462(157) MeV, and Su at 1554(65) MeV are consistent with the experimental values 
(m N = 939 MeV, m R « 1440 MeV, and m Sll « 1535 MeV). This suggests that a 
successful description of the Roper resonance depends not so much on the use of the 
dynamical quarks, but that most of the essential physics is captured by using light 
quarks to ensure the correct chiral behavior. 

To conclude, with the help of constrained curve fitting (which allows fitting of 
the excited states including the ghost rj'N states) and the overlap fermion (which 
admits calculations with realistically small quark masses), we are able to study the 
Roper and Su closer to the chiral limit than before. The physical Roper and Su are 
distinguished from the ghost two-particle intermediate states (rj'N) by checking their 
volume dependence of the weights as a function of the pion mass. We observe the 
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tendency for the Roper and Sn to cross over (with the caveat of large errors) at the 
low m n region. This is consistent with the thesis that the order reversal between the 
Roper and 5n(1535) compared to the heavy quark system is caused by the flavor-spin 
interaction between the quarks due to Goldstone boson exchanges [Hj. We support 
the notion that there is a transition from heavy quarks (where the SU(6) symmetry 
supplemented with color-spin interaction for the valence quarks gives a reasonable 
description) to light quarks (where the dynamics is dictated by chiral symmetry). It 
is suggested that this transition occurs at m n ~ 300 — 400 MeV for the nucleon. 

This work is partially supported by DOE Grants DE-FG05-84ER40154 and DE- 
FG02-95ER40907. We wish to thank A. Bernstein, G.E. Brown, L. Glozman, D.O. 
Riska, and S. Sasaki for useful discussions. 
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